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1. Introduction 

Much attention has been devoted to the spin degree of freedom in one-dimensional (ID) 
conductors, both of the linear shape (quantum wires pQ, carbon nanotubes |2j) and 
of a circular one (quantum rings Physical parameters of such systems, e.g., 

average distance between the electrons a, their total number N, their effective mass m, 
dielectric constant e, etc., can vary over a broad range or can be tuned experimentally. 
This creates unique opportunities for studying the effect of reduced dimensionality and 
strong Coulomb interactions on quantum magnetism. A number of new developments 
have generated a particular interest to the physics of a ID Wigner-crystal (WC). 
Unlike the case in higher dimensions, in ID the crossover to this strongly-correlated 
regime occurs at easily achievable electron densities (Hj, r s = a/2a,B > 4, where 
as = h 2 e/me 2 is the effective Bohr radius. Disorder has been the only major obstacle 
to realizing the ID WC experimentally [I]. A successfully solution to this problem 
has been apparently found, at least, for the case of carbon nanotubes. Very large 
r s values have been recently demonstrated in suspended nanotube devices without 
appreciable intervention of disorder effects |2] . Because of their finite length, in the 
desired range r s > 4 these devices contained only a few electrons, N < 25. Such 
finite-size systems are traditionally referred to as Wigner molecules 0. The progress 
towards realizing Wigner crystal (molecule) states in GaAs quantum wires has also 
been very encouraging P^, therefore, one may hope that they will soon follow suit. 

On the theoretical side, the ID WC is interesting because of a dramatic difference 
between the characteristic energy scales for orbital and spin dynamics. This strong 
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Figure 1. (Colour online) The instanton trajectories. (a) Schematic 
representation for N = 6 Wigner molecule on a ring, (b) The trajectories of 
1 < j < 4 electrons for N = 8 (for notations used see the main text). Inset: 
function trfl(x). The units of r and O are \/2a/so and its inverse. 

spin-charge separation has been recently predicted to cause anomalies in many 
essential electron properties, e.g., ballistic conductance [S] of quantum wires and 
persistent current of quantum rings 0. In view of the above, obtaining a reliable 
estimate of the spin-related energy scales, notably, the exchange-coupling J of the 
nearest-neighbour electrons, is desirable. It has been an outstanding challenge, though. 
As depicted in figure^a), J is determined by the acts of quantum tunnelling in which 
any two such electrons interchange. At r s 3> 1 the corresponding potential barrier 
greatly exceeds the kinetic energy of the electron pair, which makes J exponential 
small and difficult to compute numerically [7|- Attempts to derive J analytically 
(for the nontrivial case N > 2) were based on the approximation that neglects all 
degrees of freedom in the problem except the distance between the two interchanging 
electrons 8,9. We call this a Frozen Lattice Approximation (FLA) . The accuracy of 
the FLA is unclear because it is not justified by any small parameter. When a given 
pair does its exchange, it sets all other electrons in motion, too (figure To obtain 
the much needed reliable estimate of J one has to treat the spin exchange in a Wigner 
molecule (or a WC) as a truly many-body process. This is done below in this letter, 
where we compute J to the leading order in the large parameter r s . 

2. Model and results 

We assume that electrons are tightly confined in the transverse dimensions on a 
characteristic lengthscale R <C as- This allows us to treat the problem as strictly 
ID, provided we replace the Coulomb law by the appropriate effective interaction 
that goes to a finite value at distances r <§C R. We adopt a simple model form ^3 ^TJ 
U (r) = e 2 /e(r+R), which is the simplest expression that correctly captures both short- 
and long-range behaviour of the (unscreened) Coulomb interaction for any realistic 
confining potential and is similar to other forms used in the literature |12l HHj . For 
convenience, we focus on the quantum ring geometry where r = (Na/n) |sin(7ra/./Va)| 
is the chord distance and x is the coordinate along the circumference. 

In the Wigner molecule configuration electrons reside at the corners of a regular 
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Table 1. Results for Wigner molecules on a ring (finite N) and for wires (N = oo). 
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polygon. The effective low-energy Hamiltonian of such a state is given by |Zj 
21 



H=^-L 2 + J^SjSj+i+^nafa*, (1) 



3 a 

where L is the centre-of-mass angular momentum, are electron spins and n a are the 
occupation numbers of "molecular vibrations." At large r s , the total moment of inertia 
/ and the vibrational frequencies uj a are easy to compute because they are close to 
their classical values. In order to calculate J, which is more difficult, we first show that 
the asymptotically exact relation exists between J and the pair-correlation function 
(PCF) g(x) of a spin polarized ID system. For an ultrathin wire, L = ln(as/i?) 3> 1, 
it is particularly compact: 

J = (e 2 a|/2£e) g"(0), r s » 1/C. (2) 

By virtue of (J2J, the calculation of J reduces to an easier task of computing g(x). 
Using the instanton method described below we arrive at the final result 

2 

J = ; " 5/4 7— exp (-vV^Fs) , r s »l. (3) 
(2r s ) 5/4 Cea B 

The values of r\ and K are given in tabled They demonstrate that the FLA [HI El errs 
significantly in k, by about 50%, but surprisingly little in r), only by 0.7%. 

3. Three electrons on a ring 

We start with the simplest nontrivial case: N = 3. Let < Xj < 3a, j — 0, 1, 2, be 
the electron angular coordinates. We will compute the exchange coupling J between 
the j = and the j = 1 electrons. It is convenient to go to new variables: the 
relative distance of the pair, x = x\ — xq and the location of the j — 2 electron 
with respect to the centre of mass X2 = X2 — xcm — a. The motion of the centre 
of mass xcm can be ignored. We arrive at the problem with one fast (x) and one 
slow (A 2 ) degree of freedom. (Classically, X 2 = 0.) The total potential energy 
U tot (x,X 2 ) = U(x)+U [(3/2) (X 2 + a) - x/2]+U [(3/2) (X 2 + a) + x/2] has two global 
minima in the fundamental domain \x\ < 3/2a, at x = ±a, X 2 = 0. They give rise to 
the two lowest-energy multiplets: the spin-singlet ground state S + Si = with an 
orbital wavefunction $ s (x, X 2 ) and the triplet with a wavefunction $ t . Their energy 
splitting is the desired exchange coupling J . It is given by the formula (141 IT5] 

J = (2h 2 /v) JdX 2 ^d x ^\ x=0 , (4) 

where the (normalized to unity) "single- well" wavefunction &i(x, X 2 ) is the ground- 
state of the Hamiltonian with a modified potential L^ t —*U\ = U^qi (max{x, 0}, X 2 ) 
and n = to/2. Equation (0J is valid to order 0(J 2 ) with the same accuracy, the 




Figure 2. The PCF of a spin-polarized system (schematically). Regions I, II and 
III are described in the main text. 



singlet and triplet wavefunctions are symmetric and antisymmetric combinations of 
the single-well wavefunctions, <& Sj t = [$i(x, X 2 ) ± 3>i(— x,X 2 )] /\/2. 

Let us discuss the form of &i(x, X 2 ). Near its maximum at x = a, X 2 = 0, it is a 
simple Gaussian in both variables, characterized by an amplitude I of the zero-point 
motion in x and a frequency Q (a) of the zero-point oscillations in X 2 . Away from its 
maximum $1 rapidly decays at IX2I > I ^> a. This justifies the following Gaussian 
approximation! m the entire fundamental domain of x: 

$1 = (}>{x) exp [-(M/2h)n(x)X 2 2 ] , (5) 

where M = 3/i. It is important that at x <C a, where the tunnelling barrier is large, f2 
is a slow function of x. Hence, if g(x) denotes the PCF of a spin-polarized molecule, 

„ N-l 

g(x)=2 / J] dX j $ 2 t (x,X 2 ,...,X N _ 1 ), \x\<3a/2, (6) 
J i=2 

then {3} immediately entails J = (h 2 /4(i) [</>(0)/</>'(0)] g"(0). Anticipating the 
discussion below, (JSJ) is written for an arbitrary N > 2, with the notation Xj = 
Xj — xcm + (N — 1 — 2j)(a/2) being used; the PCF is normalized as appropriate in 

the WC limit, / Q 3a/2 g(x) dx = 1. 

The equations on <fi(x) and tt(x) are obtained by substituting © into the 
Schrodinger equation and neglecting terms small in I /a. This results in the dependence 
of g{x) on x sketched in figure |21 Near its x — a maximum (region III) g{x) is a 
Gaussian of width I. In the region II the quasiclassical approximation applies. Finally, 
in an ultrathin wire, C 3> 1, there is also region I, x < as, where the quasiclassical 
approximation breaks down. Fortunately, the equations on 4>{x) and f2(x) can be 
simplified there, as fl{x) ~ fi(O) and U tot {x) ~ U{x) + 2U(3a/2). Similar to |TU|. 
this leads to <j){Q)/4>' (0) ~ which, combined with the expression for J, yields 

equations J2Jl and with rj and k given by 



n = 2 



dx \ ea 
a 



1/2 



f;2 AC/tot^) (7) 



2 5 / 4 e?(0) ^/fi(a) a, 



3 

II 



^totW 



3/4 



| The Gaussian ansatz has been used previously for computing spin exchange in 3 He crystals llfi| 



(8) 
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Thus, for the N = 3 case we were able to reduce the original complicated three- 
body eigenvalue problem to routine operations of solving an ordinary differential 
equation on O(x) and taking two quadratures. The resultant rj and k are listed 
in table 2] In comparison, the FLA 8 underestimates k by about 50%. It gets rj 
correctly but only for N = 3, see more below. 

One important comment is in order. The antisymmetry of the total fermion 
wavefunction imposes certain selection rules \Tj\ f° r the allowed values of L [see Q] 
at a given total spin S. The lowest-energy L eigenstates for the two possible S values 
in the N = 3 system, 5=1/2 and S = 3/2, are \L\ = 1 and 0, respectively. Since 
J <C h 2 /I at large r s , the ground state of the system is the L = spin-quartet [71ITT|. 



4. N > 3 electrons on a ring 



In a system of more than three electrons, the single- well function <E>i(x,X) can be 
sought in the form similar to but with the argument of the exponential replaced 
by (-l/2n)(AXtM 1 / 2 )r2(x)(M 1 / 2 AX), where A/^ 1/2 = mT 1 ' 2 (1- y/2/N)/(N- 
2)]. In the language of the quantum tunnelling theory, Sl(x) is a matrix that controls 
Gaussian fluctuations AX = X — X* around the instanton trajectory X*(x), where 
X = (X2, ■ ■ ■ , Xn~i) t ■ Switching to the usual parametrization of the instanton by an 
"imaginary-time" t, we seek x(t) and X*(r) that minimize the action 



Sn 



dr 

Y 



{d T xf 



1 



(<9 T X) t M<9 T X + AC/, 



tot 



(9) 



subject to the boundary conditions x(0) = 0, x(oo) = a and X(oo) = 0. Henceforth 
L^ot is always meant to be evaluated on the instanton trajectory and AL^ t stands 
for the difference of its values at a given r and at r = 00. Repeating the steps of the 
derivation for the N = 3 case, we derive the following equations on 4>(x) and fi(x): 

{(h 2 /2fi)d 2 x - u tot (x) - (h/2)tv n(x) + e} <Mx) = o, 

where a; is a positive-definite matrix such that uS 2 — M _1 / 2 HM -1 / 2 and S is the 
matrix of the second derivatives 3jj = dx t dxj C^tot ■ The equations are mutually 
consistent if E = U^(a) + (h/2) [tr u; (a) + u>o], ^0 = fif/il 2 and fl(a) — u>(a). 

The PCF g(x) in the quasiclassical region can be written in terms of the tunnelling 
action © and the appropriate prefactor as follows: 

1/2 



9( x ) = n 



1 Q(a) huo 



5(x)-2Sjv(x) 



e(x) = / dy 



I 2 |_27rfi(x) U{x)_ 

ujq + tr SI (a) — tr Sl(y) 

. [(2/ M )AC/ tot (y)] 1 / 2 



(ii) 



(12) 



Here the action S^r is defined to be the value of the integral in © when its lower limit 
is replaced by r = t(x). For 77 we find 77 = 2Sn I ^f2r s , while k is given by equation 
JHJ after the replacement SI — ► detn. 



Letter to the Editor 



G 



5. Calculation of the instanton 

A few properties of the instanton follow from general considerations. The dimensional 
analysis of action yields Sn oc y/rl, so that 77 is indeed just a constant. Also, 
from the symmetry of the problem, Xn+i-j(t) = —Xj(t). Thus, in the special case 
of N = 3, the instanton trajectory is trivial: X2 = 0, i.e., the j — 2 electron does 
not move. This is why we were able to compute S3 in a closed form. For N > 3 
the situation is quite different: all electrons [except j — (N + l)/2 for odd N] do 
move. In order to investigate how important the motion of electrons distant from the 
j = 0, 1-pair is let us consider the N = 00 (quantum wire) case, where the far-field 
effects are the largest. If X/s were small, we could expand AL^ t in © to the second 
order in Xj to obtain the harmonic action 




where Ukuj is the Fourier transform of electron displacement %(r) = Xj — x® from 
the classical equilibrium position x® = (j — l/2)a, j £ Z, (J p (k) ~ sok ln 1 / 2 (4.15/fca) 
is the plasmon dispersion in the ID WC and sq = (e 2 /e/ia) 1 ^ 2 . Minimization of 
Sh with the specified boundary conditions yields Uj(r) oc vx®/ [{x®) 2 + v 2 t 2 ~\ , where 
v ~ (so/2) In { [(x°) 2 + SqT 2 ] /a 2 }. Substituting this formula into harmonic action 
H13fl . we find that the contributions of distant electrons to Sh rapidly decay with 
Thus, a fast convergence of 77 to its thermodynamic limit is expected as N increases. 
Encouraged by this conclusion, we undertook a direct numerical minimization of S for 
the set of N listed in table fusing standard algorithms of a popular software package 
MATLAB. The optimal trajectories that we found for the case of N — 8 are shown 
in figure ^b). As one can see, electron displacements reach some finite fractions of a 
at r = 0. This collective electron motion lowers the effective tunnelling barrier and 
causes 77 to drop below its FLA value, although only by 0.7%, see table 

Let us now discuss the prefactor k. In the inset of figure Eft)) we plot trn(x) 
computed by solving (|10|) numerically. To reduce the calculational burden, we set 
X*(t) — > instead of using the true instanton trajectory. The error in k incurred 
thereby is ~ 2%. In comparison, the FLA, where trQ(x) = const, yields k about 50% 
smaller than the correct result, similar to N = 3. 

6. Relation to current experiments 

For carbon nanotube quantum dots 2 , where the WC limit has apparently been 
realized, our formula © gives J ~ 1 K at r s — 4, which should be verifiable 
experimentally. Unfortunately, the lowest measurement temperature was 0.3 K; 
therefore, the exchange correlations may have been washed out. We hope that 
our predictions can be checked in the next round of experiments. Energy-level 
spectroscopy of quantum rings |3J |S] is another promising area where our results 
may apply. In longer ID wires, J determines the velocity v a — (tt /2)Ja/h of spin 
excitations, which can be measured by tunnelling [J, photoemission |TH], or deduced 
from the enhancement of the spin susceptibility and electron specific heat ^O]- Our 
result for v a reads (cf. tabled 

Va /v F = 5.67 (tt/£) rl^e-" 1 ^, 77 = 2.7978(2), (14) 
where vp — (n/2)(h/ma) is the Fermi velocity. 
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Note added. — After the completion of this work, we learned that Klironomos et al |19| 
independently computed t] = 2.79805(5), but not the prefactor k. These authors also 
considered a correction to i] due to a finite radius of the wire R. We can show that as 
R increases, the ratio ir/£ in (|14(l is replaced by a more complicated expression that 
tends to unity at R > as- 
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